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Abstract 

To  explain  the  increase  of  compressibility  of  vitreous 
silica  with  treasure,  it  is  assumed  that  there  is  a possibility 
of  bending  of  the  Si-O-Si  lines  in  the  structure.  Calculations 
of  network  energies  show  that  thi3  i3  a reasonable  assumption. 
Because  of  the  randomness  of  the  glass  structure,  it  is  not 
possible  at  this  time  to  get  more  than  a qualitative  agreement 
between  calculated  and  experimental  values. 

Introduction 

Most  ionic  crystals  show  a compressibility  which  decreases 
with  increasing  pressure.  This  might  be  expected,  since  in  the 
absence  of  external  pressure,  the  size  of  the  crystal  lattice 
is  such  that  the  Coulomb  forces  between  ions  which,  on  the 
whole,  act  to  pull  the  lattice  together  and  reduce  interionic 
distances  are  exactly  balanced  by  the  repulsive  forces  which 
have  appreciable  value  only  at  very  close  distances  and  whose 
magnitude  varies  inversely  as  a very  much  higher  power  of  the 
interionic  spacing  than  the  second  power  obeyed  by  the  Coulomb 
forces.  Under  increasing  pressure,  as  the  ions  are  forced 
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together,  these  repulsive  forces  increase  very  rapidly  and  it 
becomes  increasingly  difficult  to  decrease  the  spacing.  In 
other  words,  the  compressibility  decreases. 

Measurements  by  Bridgnarf  have  shown  that  the  coimresaibility 
of  silioa  glass  increases  with  increasing  pressure.  This  means 
that  it  is  unlikely  that  the  deformation  consists  entirely  of 
a decrease  of  all  interionic  distances  by  the  sam*  ratio  since 
this  would  bring  increasingly  into  rlay  the  high  power  repulsive 
forces  which  would  bring  about  a decrease  in  compressibility. 

We  must  look,  therefore,  for  some  other  mechanism.  According 
to  the  usual  ideas  on  the  structure  of  vitreous  silica  the 
structure  is  rather  open  with  numerous  empty  spaces.  This 
structure  might  allow  a straixht  line,  originally  connecting 
three  adjacent  ions,  to  bend  under  the  action  of  pressure. 

In  particular,  one  might  expect  the  line  Si-O-oi  which  according 
to  Warren’s  X-ray  work  is  or  approaches  a straight  line  to 
bend  slightly  in  the  middle.  It  was  .n  an  attempt  to  find  out 
if  this  is  a reasonable  explanation  that  the  following  calcula- 
tions were  undertaken. 

Explanation  of  Method 

While  there  is  good  evidence  for  some  degree  of  covalent 
binding  between  silicon  and  oxygen,  it  was  felt  to  be  worth- 

xRobert  B.  Sosman,  The  Properties  of  Silica,  The  Chemical 
Catalog  Co.,  Inc.,  New  York,  1927. 
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while  to  investigate  the  effect  using  a purely  ionic  model 
to  find  out  if  consideration  of  central  forcos  only  is  suffi- 
cient to  explain,  at  least  qualitatively,  the  behavior  of 
vitreous  silica  under  compression. 

As  will  be  shown,  within  the  limits  of  the  assumptions 
made,  it  is  possible  to  calculate  the  energy  of  the  network 
by  considering  the  distance  between  each  ion  and  every  other 
ion  and  the  nature  of  the  attraction  or  repulsion  between  each 
pair. 

To  find  the  equilibrium  state  at  zero  external  pressure, 
we  choose  those  values  of  the  interionic  spacings  which  make 
the  energy  a minimum.  Under  the  action  of  external  pressure, 
we  aasume  some  definite  smaller  volume  and  calculate  the  inter- 
ionic  spacing  again  to  "’.ive  the  minimum  energy  for  this  com- 
pressed state.  This  new  minimum  will  be  greater  than  the 
energy  for  zero  external  presoure;  from  the  relations  between 
this  new  energy  and  its  corresponding  volume,  it  is  possible, 
as  will  be  shown,  to  calculate  the  pressure  required  to  pro- 
duce this  state  and  hence  to  calculate  the  compressibility. 

On  the  basis  of  these  energy  calculations  alone,  it  is 
possible  to  decide,  for  each  compressed  state,  whether  minimum 
energy  is  achieved  by  maintaining  the  oi-O-Si  lines  as  straight 
lines  or  by  allowing  them  to  bend. 

Assumed  Structure 

On  the  basis  of  the  generally  accepted  structure  for  vitrocus 
silica,  it  was  assumed  that  in  the  undisturbed  state  each  Si  ion 
is  surrounded  te trahedrally  by  four  0 ions,  and 
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that  the  line  drawn  from  a Si  Ion  to  any  one  of  its  four 
neighboring  oxygens  will,  if  produced,  go  through  another 
ai  ion  which  is  itself  surrounded  tetrahedrally  by  four 
oxygens . 

On  deformation  it  i3  assumed  that  the  network  of  Si 
ions  does  not  distort  but  merely  reduces  uniformly  in  linear 
dimensions.  The  possibility  is  assumed  that  an  oxygen  ion, 
originally  lying  at  the  exact  midpoint  of  the  line  joining 
two  neighboring  Si  ions,  moves  out  of  line  as  the  Si  ions 
start  moving  together. 

Relationships  Between  Neighbors 

For  the  ourpose  of  the  energy  calculations  it  is  necessary 
i.eally  to  know  the  distance  between  each  ion  and  every  other 
:o n.  Since  in  a glassy  structure,  in  contrast  to  a crystalline 
structure,  there  is  considerable  randomness  of  structure  beyond 
the  second  or  third  nearest  neighbor,  it  is  probable  that  a 
fairly  good  a -proximation  of  tho  energy  arising  from  the  Coulomb 
interactions  can  be  obtained  b/  limiting  attention  to  a fairly 
close  groun  of  neighbors  around  each  ion,  so  chosen  that  the 
structure  outside  this  .group  is  electrically  neutral. 

In  the  undisturbed  state  the  following  neighbors  will  be 
considered,  hach  Si  ion  has  four  oxygens  at  distance  a and 
four  oi  at  distance  b where  b is  equal  to  2a.  Further  out 
there  will  be  twelve  oxygens  at  a distance  d where 
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Of  these  we  consider  only  six  which  means  that  the  net 
charge  on  the  neighbors  is  *4  which  exactly  balances  the 
charge  of  plus  4 on  the  silicon  and  leaves  a neutral  mass 
surrounding  these  neighbors.  The  contribution  of  this  neutral 
mass  is  assumed  to  be  zero  on  the  average. 

Around  each  oxygen  there  are  two  closest  silicon  ions 
at  a distance  a and  six  oxygens  at  a distance  c where 
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Of  these  silicons  at  distance  d we  count  only  3/2  in 
order  to  maintain  neutrality  outside  these  neighbors.  On 
distortion  we  do  not  assume  that  b is  equal  to  2a  but  do 
assume  that  on  the  average 


2 n 6 
c = ~3 — a 


(3) 


This  is  equivalent  to  saying  that  after  distortion  each 
silicon  is  still  surrounded  tetrahedraily  by  four  oxygens 
although  the  Si-O-bi  lines  are  no  longer  straight. 

The  quantity  d will  now  involve  both  a and  b.  On  the 


average  it  will  be  assumed  that  d is  given  by 

d8  = aa  + ba  ♦ 2ab 


(4) 


This  averaging  of  distances  to  second  nearest  neighbors 
seems  justified  because  of  the  randomness  of  a glassy  structure 
as  contrasted  to  the  definite  periodicity  of  a crystalline 
structure , 
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Enerca 

The  energy  of  the  lattice  arising  from  the  Coulomb  inter- 
actions is 


\ 1 Z*>  (5) 

where  e is  the  magnitude  of  the  charge  on  an  electron,  Zn  is 

the  number  of  plus  charges  carried  by  ion  n and  Cf^is  the 

potential  at  the  location  of  ion  n arising  from  all  its 

neighbors.  The  1/2  arises  becauses  in  this  summation  every 

ion  oair  is  considered  twice. 

The  energy  per  Si08  molecule  will  therefore  be 

\ fu»  + 2(-2o)  (6) 

where  (f  is  the  average  value  of  the  potential  at  a silicon 

I 

location.  Thi3  nay  be  written 


2e(:(s-  <p.  ) (7) 

On  the  basis  of  our  previous  assumptions  concerning  neighbors 

(jSv  = -4  2e  + 4 ^ . 6 2e 

b "d” 


« e,  8 + 16  - 12  x 

V — T~  ) 


= 2 i£0  - 6 2e^  + 2.  ^ 
° c 2 u 


a 


~ Q 1 8 - 12  + 6 x 

(a  c“  7 } 


(8) 


This  gives  for  the  energy  per  Si08  molecule 

2e®  rl6  +16+12  - 18  x 

"F  c"  T) 

o®  (-32  + 22  + Q.  - J6  J 

a b c d 


or 


(9) 
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Contribution  of  Repulsive  Terms 

The  above  magnitude  becomes  increasingly  negative  as 
the  network  shrinks  in  size.  A minimum  value  for  the  energy 
is  possible  because  of  the  repulsion  between  ions  which  arises 
at  small  distances  of  separation  and  whose  magnitude  is 
frequently  approximated  by  an  expression  of  the  type 
A 

r*'  (10) 

whererlmay  be  in  the  neighborhood  of  10.  Since  this  falls 
off  so  rapidly  with  distance,  not  too  much  error  ^s  introduced 
in  considering  interactions  only  between  nearest  neighbors. 

In  this  case  the  energy  per  molecule  arising  from  this  source 
is  of  the  form 


e*B 


where  B is  a constant  and  a is  the  silicon-oxygen  distance. 

The  factor  e*  is  included  in  the  constant  merely  for  convenience 
in  combing  this  repulsive  energy  with  the  Coulomb  energy  given 
in  equation  (9).  The  repulsive  energy  (11)  added  to  the 
Coulomb  onorgy  (9)  has  a minimum  value  for  some  valuo  of  a 
which  is  tho  zero  pressure  equilibrium  value  of  a. 

Relations  Bctv/ecn  Enemy  and  Prcssuro 

As  shown,  for  instance,  in  Chapter  II  of  Seitz'  Modern 
Theory  of  Solids8  tho  relation  between  the  pressure,  energy 
and  volume  of  a solid  is 


2 Frederick  Seitz,  The  Modern  Theory  of  bolids,  McGraw-Hill  Book 
Co.,  Now  York,  1940 
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Under  zero  pressure  the  solid  reaches  a state  of  minimum 
energy  where 

4=° 

The  connection  between  energy  and  volume  as  deduced  from 
our  previous  calculation  enables  us  to  assign  a pressure  to 
each  volume  state  of  the  glass. 

To  compute  the  compressibility  in  any  state  we  uco  the 
relation 

where/?  is  the  compressibility.  In  this  formula  E is  the 
energy  per  mole  and  V is  the  molar  volume.  The  connection  be- 
tween V and  the  quantity  b we  have  been  using  is  obtained  by 
using  the  fact  that  V is  proportional  to  the  cube  of  b and  the 
fact  that  we  know  from,  experiment  the  molar  vol’ure  of  vitreous 
silica  m the  unstressed  state. 

Factors  Determining  the  Form  of  N etwork 

In  the  following  calculations  it  wilo.  be  assumed  that  as 

the  network  is  co  pressed  tnere  1 s no  distortion  if  the  part 

of  the  network  which  consists  only  of  oi  ions  but  merely  a 

shrinkage.  There  is  allowance,  however,  for  a bending  of  the 

^i-O-Si  line3  and  this  is  handled  as  follows. 

A value  of  the  quantity  b is  assumed  which  is  the  closest 

Si-Si  spacing.  Values  of  D anc.  n for  tli  constants  in  the 

repulsion  energy  term  are  also  assumed.  Then  for  each  value  of 

b a 3eris3  of  values  of  a is  taken  whose  smallest  value  is  b 

2 
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and  the  energy  of  the  network  is  computed.  That  value  of  a 
which  makes  the  energy  a minimum  for  the  value  of  b chosen  is 
taken  as  the  correct  one.  Thus,  the  bending  or  non-bending 
of  the  Si-O-Si  lines  (b  less  than  2a  or  b equals  2a)  is  merely 
a matter  of  what  configuration  gives  the  lowest  energy  for  the 
value  of  b selected. 

Results 

As  a start  it  was  assumed  that  the  Si-Si  spacing  In  the 
unstressed  state  is  3.21+.  A.  For  want  of  any  very  good  informa- 
tion the  assumption  was  made  that  silicon  and  oxygen  ions  repel 
each  other  according  to  an  inverse  eleventh  power  law  or  that 
the  energy  arising  from  this  effect  varies  as  the  inverse  tenth 
power  of  the  separation. 

First  it  was  assumed  that  3 has  that  value  which  makes 
the  energy  a minimum  for  a silicon-silicon  distance  of  3,2k  A. 

W hen  this  value  was  used,  it  was  found  that  for  the  entire 
range  of  compressions  covered  in  Bridgman’s  work,  the  energy 
relations  required  that  the  oxygens  lie  along  the  silicon- 
silicon  line.  This,  as  would  be  exoected,  gave  a coimores slbility 
which  slowly  decreased  with  increasing  pressure.  The  Si-O-Si 
lines  did  not  bend  until  the  volume  decreased  considerably 
greater  than  that  actually  investigated.  In  units  such  that 

eaB 

10 

a 

has  the  dimensions  of  ergs  ner  Si08  molecule  this  value  of  B is 
about  120, 

The  explanation  of  this  is  very  probably  that,  because  of 
the  somewhat  random  manner  in  which  vitreous  silica  or  any  gla~o 
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ls  put  together,  the  Si-O-Si  lines  are  not  straight  to  begin 
with  as  they  would  be  in  the  high  cristobalito  structure  but 
are  already  slightly  bent  and  ready  to  bend  ■"ore  immediately 
when  pressure  is  applied.  This  means  that  the  unstressed 
stato  does  not  correspond  to  a minimum  energy  state  within  the 
framework  of  the  simple  assumptions  made.  The  effect  of  the 
random  structure  of  the  gla^s  j.o  , in  fact,  el-'ost  as  j.f  the 
ideal  strvctu.ro,  at  first  assumed,  were  subjected  even  at  zero 
presov.-c  t'  a proscure . sufficient  to  bend  the  oi-O-Si  linos 
slightly. 

The  network  was  then  treated  as  _f,  at  tho  equilibrium 
spacing  of  3.21;  A,  it  w«re  unber  -.uch  a pressure  that,  on 
increasing  tho  -res sure,  the  Si-O-Si  lines  would  start  bending 
within  the  range  covered  by  the  experiments. 

A value  of  160  for  B required  that  even  at  a oi-Si  spacing 
of  3.21;  A tho  Si-O-Si  lines  we~e  already  bent.  A value  of 
lA  meant  that  the  Si-O-Si  lines  remained  straight  on  compression 
until  the  Si-oi  distance  was  reduced  from  3.21;  A to  3.221;  A; 
on  further  co -press ion,  minimum,  energy  was  achieved  by  displacing 
the  oxygon  ion  from  thi3  line  but  keening  it  An  tuo  plane  which 
perpendicularly  bisected  the  lino. 

i?rom  3.21;  until  the  line  starts  to  bond,  there  is  a 
slight  decrease  m compressibility.  At  the  onset  of  bending 
there  is  a s^ddon  discontinuous  change  in  compressibility 
to  a considerably  higher  value  which  on  further  compression 
slowly  starts  to  oecreaso.  Th~  suddonness  of  the  increase 
in  compressibility  is  not  in  agreement  with  experiment.  Tho 
gradualness  actually  found  is  again  probably  caused  by 
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the  randomheafc  of  the  vitreous  structure.  This  randomness 
probably  results  in  a different  dogree  of  bending  of  the 
various  Si-O-Si  lines  throughout  the  structure.  Increase 
in  compressibility  arises  only  when  tho  linos  start  bonding. 

It  would  then  seem  reasonable  to  assume  that  in  tho  normal 
state  some  of  these  lines  are  straight  and  some  slightly  bent. 
Increase  in  pressure  increasos  the  number  actually  bent,  and 
this  gradual  onset  of  bonding  gives  rise  to  tho  slow  rather 
than  sudden  increase  in  compressibility. 

Conclusions 

Since  compressibility  of  vitreous  silica  starts  to 
increase  immediately  on  increase  of  pressure,  the  calculations 
would  indicate  that  under  aero  or  atmospheric  prossure 
there  is  already  some  bending  of  the  Si-O-Si  lines. 

Because  the  increase  in  compressibility  is  gradual 
rather  than  suddon,  it  might  be  concluded  that  there  is 
some  randomness  in  tho  degroe  of  bonding  of  tho  Si-o-Si  linos. 

Those  conclusions  aro  basod  on  calculations  using  a 
purely  ionic  model,  and  no  attempt  has  been  made  to  ostimate 
what  the  effect  of  some  degree  of  covaloncy  in  the  Si-0 
binding  might  be. 
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